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ABSTRACT 

We report on submillimetre bolometer observations of the isolated neutron star 
RXJ1856. 5-3754 using the LABOCA bolometer array on the Atacama Pathfinder 
Experiment (APEX) Telescope. No cold dust continuum emission peak at the po- 
sition of RXJ1856. 5-3754 was detected. The 3a flux density upper limit of 5mJy 
translates into a cold dust mass limit of a few earth masses. We use the new sub- 
millimetre limit, together with a previously obtained H-band limit, to constrain the 
presence of a gaseous, circumpulsar disc. Adopting a simple irradiated-disc model, we 
obtain a mass accretion limit of M <i 10^^ g s~^, and a maximum outer disc radius of 
~ 10^^ cm. By examining the projected proper motion of RX J1856. 5-3754, we specu- 
late about a possible encounter of the neutron star with a dense fragment of the CrA 
molecular cloud a few thousand years ago. 
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1 INTRODUCTION 

Since the discov ery of planets a round the pulsar 
PSR 1257+12 by IWolszczan fc Fraill (|l992D . dusty discs 
around pulsars have become interesting to observers, who 
have been trying to detect them in the infrared (IR) or at 
(sub-)millim eter wavelengt h s. Mo st co mprehensive are the 
search es by iLohmer et al.l l|2004h and iGreaves fc Holland! 
(|2000l ). Previous searches had concentrated mainly on 
recycled, thus formerly accreting, radio pulsars as the most 
likely objects to be surrounded by dusty discs. Neither 
of the above mentioned surveys detected dust emission 
around the planet-hosting PSR 12574-12, nor did the 
recent search with Spitzer at 24 fim and at 70 /im by 
iBrvden et al.l (|2006l ). However, dusty discs could also be 
present around non-recycled neutron stars (NSs). Following 
the supernova explosion, which creates the NS, some of 
the explosion ejecta may fail to escape and remain bound 
- forming a fallback disc. Such fallback d iscs are a general 
prediction of current supernova models jMichel fc Desslen 
Il98ll : [Chevalieij|l989l ). and have been invoked by a number 
of author s to explain a variety of phenomena related to 
NSs (e.g.. IChatteriee et al.|[200ol : lAlpaj|200ll : iMenou et all 
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I2OOII: iBlackman fc Pernal 120041 ). Recently, IWang et al] 
(2006|) reported the discovery of mid-infrared emission 
from a cool disc a round the An omalous X-ray pulsar AXP 
4U 01424-61. Wang et al.l bOOdt) inte r preted their detection 
as a passive disc, while lErtan et al.l (|2007| ) argued that it 
could also originate from an actively accreting disc. To 
date, AXP 4U 0142-1-61 is the only isolated neutron star 
for which a fallback disc is believed to have been detected. 
Such discs appear to be rare. According to lEksI fc Alparl 
(|2005h fallback disc are rare because they are likely to be 
disrupted when the newly born NS spins rapidly through 
the propeller stage, at which in-flowing matter, instead of 
being accreted, would be expelled. The fallback discs can 
survive if the initial NS spin is slow e nough (^ 40 ms at a 
magnetic moment of ^ = 10^'^ G cm^). |jone3 (|2007l ) studied 
the effect of pulsar wind induced ablation of fallback discs. 
He concluded that long-lived discs could be present in 
many pulsars without exceeding published limits on IR 
luminosity. 

In the following, we report on submillimetre observations 
of RX J1856. 5-3754, which is the brightest and closest 
member of a class of NSs neglected so far in the search 
for circumstellar discs, the X-ray thermal isolated neutron 
stars. They are peculiar because they show pure thermal 
soft X-ray spectra without any (confirmed) non-thermal 
emission, especially no confirmed radio emission; for re- 
views, see iKaplaui (|2008. ) or .Haberf (2007. ). These objects 
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have periods in the range of 4 to 12 s. Thus, they are 
much slower than the bulk of radio pulsars. Their X-ray 
pulse periods and period derivatives are similar to those of 
the AXPs and soft gamma-ra y repeaters, SGRs, (see, e.g., 
iKaplan fc van Kerkwii^ 120091 '). This has led to discussions 
about wh ether t h ey ni a y be related to those objects. For 
example, lAlpaJ (|200ll . 120071 ') suggested that the X-ray 
thermal isolated NSs may simply have accretion discs with 
smaller masses than those of the AXPs. We note that it is 
currently not clear whether AXPs host accretion discs. The 
prevailing model for AXPs and SGRs is the magnetar model 
- isolated, young neutron stars wi th exceptionally high 
(«10"G) surface magnetic fields jPuncan fc Thompson! 
Il992l ). 

RX J1856. 5-3754 has been observed at radio wave- 
lengths (to date without detection), near infrared (without 
detection), optical and X-ray wavelengths. The striking 
absence of deviations from a pure blackbody in its X-ray 
spectrum has led to lively discussions on the nature of 
the object and of its X-ray emission, including highly 

IHo et al.i [2007). condensed 



magnetized at mospheres (e 
surfaces (e.g . 



iDrake et al 



.Burwitz et al.l 120031 ') nd quark stars (e.g.. 
20021). Th e most recent, preliminary, par- 



allax by Kaplan et al.l l|2007l ) is 1671 ijg pc and marks 



RX J 1856. 5-3754 as o n e of t he closest NSs currently known. 
iTiengo fc MeregheTtH (|2007l ') discovered the ~ 7 s pulsa- 
tions of RX J1856. 5-3754, which has an extremely small 
pulsed fraction in the 0.15-1.2 keV range o f only 1.2%. A 
timing study by Ivan Kerkwiik fc KaplanI l|200^ inferred 
a magnetic field of 1.5 x 10^^ G assuming spin-down by 
dipole radiation. Interestingly, the spin-down luminosity 
is not high enough to explain the Hq nebula around 
RX J1 856. 5-3754, discovered by Ivan Kerkwiik fc Kulkarriil 
l|200ll ). The nature of the Ha nebula remains unclear 
to date (see, e.g., Ivan Kerkwiik fc KaplanI |2008| ). These 
authors also derived a characteristic age of 4Myrs, which 
is much older than the kinematic age, 0.5 Myrs, obtained 
assu ming an origin in the Up per Scorpius OB association 
(e.g.. IWalter fc Lattimei||2002l '). At such ages, any fallback 
disc around the source would be expected to be rather cool. 



"Wi thin 
I'Walter 



dista nce 

LattimeJ (|2002l ') 



117 ± 12 PC by 



error bars, 

to 167litg pc by iKaplan et al.l 

I 2007h . RXJ1856. 5-3754 appears to be in the outskirts of 
the Corona Australis star-forming cloud, whose distance 
is relatively well known: 129 ± llpc from the orbit so- 
lut ion of the double- lined spectroscopic binary TY CrA 
by ICasev et al. (Il998f): for a detailed discu ssion we refer 
to Sect. 1. 3 in iNeuhiiuser fc ForbrichI 1I2OO8I '). Bondi-Hoyle 
accretion (|Bondi fc Hovldll944l : lBondilll952l ') from the inter- 
stellar medium (ISM) is unlikely to be a major contributor 
to the observed X-ray luminosi ty today, given the high 
velocity of this neutron star (|van Kerkwiik fc Kulkarnil 
I2OOII ). Furthermore, it is likely that isolated NSs accrete 
at sub-Bondi rates (|Perna et al.ll2003 ). The accretion rate 
scales inve rsely with the magne ti c moment as At ~ ^'^ a-c - 
cording to iToropina et all (|2006l ) ; iRomanova et al] (|2003l ). 
RXJ1856. 5-3754 has a high m agnetic field, 1.5 x lO" G 
(|van Kerkwiik fc KaplanI l2008l l. reduci ng the possibility 
of acc retion even further. However, as iDrake fc Marshall 
l|2003l ) noted, RX J1856. 5-3754 might have passed more 
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Figure 1. The IRAS/IRIS 100/im band of the neighbourhood 
of RX J 1856. 5-3754. Strong infrared emission, signaling the pres- 
ence of dust, is shown by darker areas. The peak dust emission 
outlines the CrA star forming region. The position of RX J1856.5- 
3754, as well the direction of its proper motion, is indicated. The 
dashed white line shows the rough projection of the proper mo- 
tion into the past. About 3000 years ago the neutron star might 
have crossed a region of the denser dust filament - depending on 
the actual distances of the filament and RX J1856. 5-3754. The 
boxes 1 and 2 cover the same area, see section [4.3l for discussion. 



dense regions than it now resides in. Indeed the IRAS 
100 ^m image shows some cloud fragments projected along 
the past trajectory of RXJ1856. 5-3754 (see Fig.[T}, opening 
the possibility that the neutron star collected material in 
the past. 

Here, we investigate the surroundings of RX J1856. 5-3754 
for cold dust as one might expect either in an old fallback 
disc or in a dense gas disc. The gas content of a fallback disc 
is not well known. The formation from supernova ejecta 
would be expect ed to produce a disc m ainly consisting of 
heavy elements. ICurrie fc HansenI (|2007| ) showed that gas 
from viscous, circumpulsar discs typically depletes on very 
short, « 10^ yr, timescales due to high temperatures during 
the early evolutionary stages, subsequent rapid spread of 
the disc and cooling allowing fo r the heavy elements t o con - 
dense onto grains. In addition, iPhillips fc Chandled (ll994D 
have shown that for a pulsar moving with « lOOkms"^ 
through the ISM, small dust grains with radii of ^ 0.1 fim 
are spiraled in from a circumpulsar disc to the neutron star 
on a time scale of order lO'^ years. Thus, a circumpulsar 
disc at the age of RX J1856. 5-3754 is likely to be similar 
to debris discs, which a re dust-rich and d ominated by 
larger grains (see, e.g., iKrivov et ah! |2009| ) as opposed 
to protostellar discs, which are gas-rich. If, on the other 
hand, RX J1856. 5-3754 has collected new, gas-rich material 
from the outskirts of the Corona Australis star-forming 
cloud, the disc composition could more closely resemble a 
protostellar disc. Given the uncertainties in the gas content 
and composition of a possible disc around RX J1856. 5-3754, 
we will consider here both types of discs - dense, gas-rich 
protostellar-type discs, and gas-poor, dust-rich debris-type 
discs. 



Submillimetre observations of RX J1856. 5-3754 3 



2 OBSERVATIONS AND DATA REDUCTION 

RXJ1856. 5-3754 was observed from the 18th to the 24th 
of April 2008 using th e 295-element Larg e Apex Bolome- 
ter Camera (LABOCA. ISiringo et al.]|2009l ) on the Atac ama 
Path finder Experiment (APEX) telescope (jciistcn e t al.l 
l2006l 'l. LABOCA operates at 345GHz (870^im) with a band- 
width of 60 GHz. The total on-source observing time was 
18 hours. Mapping was performed on-the-fly using spiral 
patterns to ensure even coverage of the source. To obtain 
the lowest possible flux limit, the most sensitive part of 
the array was centered on the source. The observing con- 
ditions were good to excellent with precipitable water va- 
por levels typically below 0.5 mm. Sky-dips were performed 
hourly, and combined with radiometer re adings to ob t ain ac - 
curate opacity estimates as described bv lWeifi et al.l (|2008l ). 
Th e absolute flu x calib ration follows the method outlined 
by ISiringo et al.l (|2009l ), and is expected to be accurate 
to within « 10% (e.g., ICreve et all 120091 '). The nominal 
LABOCA beam is FWHM 18.6"±1.0", and the pointing un- 
certainty is ~ 4''0. The data were reduced using the Febru- 
ary 2008 release of the BoA reduction packagqj (SchuUer et 
al. in prep.). Reduction steps included the removal of corre- 
lated (atmospheric and electronic) noise in several iterations, 
and the removal of spikes and excessively noisy channels. 

In the second part of the last observation night, we 
found a sharp drop in the root mean square (rms) of individ- 
ual bolometer time streams. As these time streams did not 
contain the expected correlated sky component, these obser- 
vations (about two hours) were excluded from further con- 
sideration. To filter out extended emission and residual at- 
mospheric noise at low frequencies, we used direct high-pass 
filtering of the raw time streams applying the fast Fourier 
transform (FFT). This was done by rejecting frequencies 
below a chosen cutoff, uc- 

For each scan, a map with 3" x 3" sized pixels was 
constructed, weighting the data by the inverse rms of each 
reduced time stream. All maps were then co-added, and 
the final map was smoothed with the LABOCA beam, 
18.6". The map rms per beam was computed as the pixel 
rms in the smoothed map. Fig. [2l applying the BoA-task 
computeRms . 



3 RESULTS 

Fig. [2] is a filtered image, where all modes corresponding to 
spatial sizes larger than about 150" are filtered out (the 
median scan speed is around 96 arcsec s~^). The same 
area as visible in Fig. [2] was considered when we com- 
puted the rms/beam for filtered or unfiltered maps. No 
submillimetre-emission has been found at or near the po- 
sition of RX J1856. 5-3754. The rms has been measured as 
3.02 mjy/beam for the map obtained without direct FFT 
filter, and 1.49 mJy/beam for the direct-FFT-filtered map. 
Applying the filter on a test point source of 10 mjy results 
in recovering of 89 % of the original flux. We correct for this 



^ More on LABOCA calibration can be found at: 
| http: / /www. apex- telescope.org/bolomcter/laboca/calibration/ 
^ httpT// www. apex- telescope.org/bolometer/laboca/boa/ 
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Figure 2. Tliis LABOCA-345 GHz-continuum-map has been fil- 
tered for all modes corresponding to spatial sizes larger than 
about 150". Darker areas indicate higher submillimetre flux. The 
circle is centered on the position of RX J1856. 5-3754 with a ra- 
dius of 18.6" (double beam size). No obvious submillimetre point 
source has been detected at the position of RX J1856. 5-3754. This 
map has a rms of 1.49 mJy/beam, corresponding to a 3cr flux den- 
sity limit of 5 mJy (see text). 



factor and obtain an 3a limit of 5mJy/beam for the filtered 
map of RXJ1856. 5-3754. 



4 DISCUSSION 
4.1 Dust mass limit 

In contrast to infrared emission, submillimetre emission is 
nearly always optically thin and traces the total dust mass. 
The conventional formula for calculating this dust mass is: 



(1) 



where u is the frequency, is the flux density, D is the 
distance, is the Planck function at a dust temperature 
Td, and k'^^ is the dust mass absorption coefficient. 

The value of k'^i^ depends on the grain parameters. 
As noted in the introduction, we consider in the following 
dust in protostellar discs as well as in debris discs to re- 
flect the wide range of possible disc compositions around 
RX J1856. 5-3754. Since the gas content in fallback discs is 
unconstrained, we calculate the dust mass limits for both 
disc types. The overall mass of a protostellar-composition 
disc would be larger by roughly a factor 100, the conven- 
tional gas-to-dust-ratio. 

In the case of protostellar-disc dust, small grains are 
present. For such protostellar dust , we use kI^v — 3 c m^ g~ ^ 
according to the relation by iBeckwith" et al.l (|l990l ): 
Ki, = 0.1(!//10^^ Hz)'' cm^ g"'^ with 1 3 = 1, and where Ki, is 
the opacity of dust and gas together. IBeckwith et al.l (ll99Gf) 
assumed a gas-to-dust-ratio of 100. Ossenkopf fc Henningj 
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Table 1. Dust mass limits for the measured 870 fim 3cr flux den- 
sity limit of 5mjy at different dust temperatures, Td and mass 
absorption coefflcients, h'^345GHz- 



dust type 






Mjj -limit 




[ cm^g-i ] 


[K] 


[M®] 


Protostellar 


3 


17 


2.0 


Debris disc 


0.3 


30 


9.1 


Debris disc 


1.7 


17 


3.6 


Debris disc 


1.7 


30 


1.6 


Debris disc 


1.7 


100 


0.4 



l| 19941 ) estimated uncertainties for ice-covered dust srains 
in protostellar molecular cloud cores to be usually within 
a factor of 5. In the case of debris discs, reported dust 
mass absorption coefficients at submilli metre wave l ength 
range from w'^^ = 0.3 - l.Tcm^g'^ (jPent et al.l I2OO0I : 
iNaiita fc Williams! |2005| ). 

Dust temperatures, Td, can range from « lOK for cold 
interstellar dust to ~ 1500 K (sublimation temperature of 
silicate dust). Due to missing constraints in the FIR we as- 
sume as a first approximation that the bulk of the dust visi- 
ble at 870/xm emits at ar ound 17 K. A c cordin g to the X-ray 
irradiated disc model bv lVrtilek et"al] (|l990D . and with the 
X-ray luminosity of RX J1856. 5-3754, Lx = 8.8 x 10^^ erg 
s~^ at a distance of d = 167 pc jPrake et al.ll200^ ). this tem- 
perature is reached at a radius of « 7 x 10^'^ cm from the 
neutron star. This value is within the disc sizes di s cusse d 
in models of NS fallback discs, e.g., bv lErtan etall l|2007l ). 
Dust in debris discs is thought to be warmer than the one in 
protostellar discs, with a range fr om 10 K up t o a few hun- 
dred K, but mostly below 110 K (|Wvattll20oi ). In Table [1] 
we give dust mass limits for the different sets of parameters. 

Another uncertainty of the mass estimation is related 
to the porosity of the dust. Porous grains are expected 
both in the cold ISM as well as in debris discs. Their in- 
fluence on k'*^ and Td can result in masses of a factor three 
smaller with respect to mas ses derived for non-porous grains 
l|Voshchinnikov et al.|[2006l '). In case of fast-moving neutron 
stars like RX J1856, 'sandblasting' by particles from the ISM 
fragments the dust (IPhillips fc Chandlerlll993) and reduces 
the number of porous grains. 

4.2 Constraints on the presence of a gaseous disc 
component 

If discs around isolated NSs are (still) gas-rich, they are 
expected to be emitting mostly in the IR and longer 
wavelengths . In t he following we use the disc model by 
iPerna et all (I2OOOI ) to constrain the presence and properties 
of a potential gaseous circumstellar disc around RXJ1856. 
The model is based on an active, viscous disc, and takes irra- 
dia tion through X-rays into account following the approach 
by IVrtilek etall (|l990l ). We note that the model does not 
include radiative transfer modeling, and it assumes that the 
disc is optically thick throughout. However, given our scarce 
knowledge about neutron star discs in general, and the lack 
of an actual IR or mm detection for RX J1856 we restrict to 
this model in the present work. 

Given the very low X-ray luminosity of RX J1856 (8.8 x 



10^^ erg s"^ at 167 pc, iDrake et all l2002h . if there is an 
active disc, it must be in the propeller phase, which is 
equivalent to saying that the magnetospheric radius, Rm ~ 
6.67 X Af /(lO^^g s~^) cm, must be smaller than the coro- 
tation radius, Rco ~ 1.5 x 10**(P/s)^''^ cm, where P is the 
period of the X-ray pulsar. Therefore, in our modelling of the 
disc emission, the inner radius of the disc, set equal to Rm, 
is restricted to be smaller than Rco- With the mag netic field 
streng th estimate 1.5 x 10^"^ G by Ivan Kerkwiik fc KaplanI 
l|2008h . this yields the condition that only accretion rate 
values M < 10^^ g/s are allowed (or else the star would be 
accreting and hence be much brighter in X-rays). 

In addition to our submillimetre observations, we use 
previously acquired data at near-infrared and infrared wave- 
lengths. More specifically, we use the H-hand limit derived 
at the position of RX J 1856, H = 21.54 ± 24 mag, cor- 
respo nding to a fiux density of 0.0025 mjy jPosselt et al.l 
I2OO9I ). Furthermore, we derive the conservative IRAS /IRIS 
fiux density limits at th e position of the neutron star apply- 
ing the IRIS atlas maps (|Miville-Deschenes fc Lagachellioosi : 
iNeugebauer et al.l ITisl ). These limits,< 1.9 Jy, < 2.7 Jy, 
< 2.4 Jy, and < 9.6 Jy at 12, 25, 60 and 100 microns re- 
spectively, are too large to constrain our model. 

Within the allowed range of M, we then investigated 
the expected disc emission as a function of M, and the disc 
boundaries Ri„ and Rout- We find that Rin infiuences the 
emission at wavelengths shorter than the ones considered 
here, and therefore the model parameters that can be con- 
strained are M and Rout (assuming a typical disc inclination 
of 60°). We find that the H-band limit restricts the accre- 
tion rate to be M ^ 10^* g/s, largely independent of Rout- 
The submillimetre Scr limit, on the other hand, restricts the 
outer radius to be Rout ^ lO^** cm or 7 AU. This outer ra- 
dius is essentially independent of the accretion rate since the 
emission in the submillimetre is dominated by X-ray irradi- 
ation. Fig. [3] shows the predicted disc emission for a range 
of outer disc radii, and the maximum allowed value of the 
accretion rate, M 10^'* g/s. For the kinematic age of the 
pulsar, 0.5 Myrs, and the derived mass accretion limit, we 
can estimate the accretion rate AI{to) during the initial time 
to ~ 1.5 X lO"'^ Rdj{to) yr (Menou et al. 2001; Rd,7 is the 
initial disc radius in units of 10^ cm) before the power law 
decay begins (Cannizzo, Lee fc Goodman 1990), and hence 
the initial disc mass M{to) - M(io)io 7 X 10'^ Mq. This 
upper limit for the initial disc mass is only a very small frac- 
tion of the overall fallback mass discus s ed to be in the range 
of ~ O.OOIMq to O.IMq bv lLin et al.1 (|l99ll ) and lChevalie"il 
(i 198911 . 



4.3 Looking at the past 

As noted above and shown in Fig. [TJ the projected trajec- 
tory of RX J1856. 5-3754 crossed a CrA molecular cloud frag- 
ment roughly 3000 years ago. The current distance estimates 
for the CrA cloud 130 pc: iNeuhauser fc ForbrichI l2008h 
and RXJ185 6 . 5-375 4 (167ljg pc; preli minary distance by 
iKaplan et al] (120071 ) or 117 ± 12 pc by IWalter fc Lattimed 
(200^) are close enough to each other to speculate about a 
possible e ncounter bet w een t he neutron star and the cloud 
fragment. IChini et al.l (|2003D investigated the R Corona 
Australis molecular cloud in detail by submillimetre contin- 
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Figure 3. The predicted flux density for an optically thick, 
gaseous disc around RXJ1856. The mass accretion is set to 
M = 10^* g s~^ (maximum allowed value by the IR limit), and 
the different curves show the emission for a range of outer disc 
radii allowed by the sub-mm limit. Plotted are also the flux den- 
sity limits from our APEX submillimetre observations and the 
i?-band. 

uum and line observations. Unfortunately, their map does 
not cover the cloud fragment west of RXJ1856. For the 
cloud cores that they investigated in the extended eastern re- 
gion, they derived densities on the order of 5 x 10~^^ gcm""^ 
to 3 X 10""gcm"^. These are a factor 10^ to 10'^ larger 
than what is usually assumed for the ISM. If RXJ1856 
crossed such a high density patch, this could have af- 
fec ted first a possible surrou nding fallback disc, destroying 
it (jPhillips fc Chandleil[l99l ) and second, it could have en- 
abled the neutron star to accrete at higher rates from the 
ISM. The extended emission seen in the IRAS image un- 
fortunately only allows us a crude estimation of a lower 
density limit in the fragment: We consider the IRAS/IRIS 
100 nm emission in box 1 centered at the fragment in com- 
parison to the emission in a box 2 of the same size, but 
centered around an area devoid of enhanced IRAS emission 
(see Fig. [T]), that supposedly represents an average density 
value for the diffuse interstellar medium. By working with 
the flux density ratio we can ignore the problems involved in 
de-convolving the extended emission from the instrumental 
response. The flux density ratio, translating into a column 
density ratio, is 2.6. The distance to the cloud/fragment is 
at least 130 pc and we take this as the minimum distance, 
d2 , along which the column density, = n2 x d2 , is derived 
from the number density, n2, in box 2. The angular extent 
of the fragment along the direction of the proper motion of 
RXJ1856 is around 18', corresponding to ~ Ipc at a dis- 
tance of 130 pc. We assume the same extent of the fragment 
along the line of sight, di ~ 1 pc; the column density, A^i, 
in box 1 is A'^i = ni X di -I- n2 X d2. Therefore, the ratio of 
the number density in the fragment to the number density 
of the diffuse ISM is at least 200, ni/n2 ^ 200. Accord- 



ing to iPhilhps fc Chandler] | 19941 ). at a number density of 
200 particles cm~^, the smallest dust grains (< 0.1 /^m ) in 
a fallback disc can be removed within 3600 years, assum- 
ing a neutron star velocity of 100kms~^ as e stimated for 
RX J1856 (e.g.. Ivan Kerkwiik fc KulkarnillioOll ') . Larger disc 
particles, however, would still survive and can be detected in 
principle in the submillimetre wavelength regime. The fact 
that we do not detect them in our APEX observation could 
mean that RXJ1856 passed a patch of much higher ISM 
density, destroying its potential fallback disc. However, since 
the data about a possible encounter between the cloud frag- 
ment and RXJ1856 is currently insufflcient (unknown 3D 
motion of the neutron star and distance of the fragment), 
it is also possible that there never was a circumstellar disc 
with a dust mass higher than few earth masses. 

Overall, the above speculations show that a better 
knowledge of the past neighbourhood of RX J1856 is desir- 
able. This can be obtained in principle by further millimeter 
continuum or molecular line observations. 



5 CONCLUSIONS 

We investigated RX J1856. 5-3754 for 870 /im continuum 
emission which would be indicative of a cold dusty disc 
around the neutron star. The derived deep flux density 
limit translates into a dust mass limit of few earth masses. 
Applying the irr adiated (gas-rich) accretion disc model by 
iPerna et al.l l|200d ) , together with further observational con- 
straints, we obtained a mass accretion limit of M ^ 10^* 
g s~^, and a constraint on the outer disc radius, Rout, to 
be smaller than 10^* cm or 7AU. Looking at the projected 
proper motion of RX J1856. 5-3754, we note that the neu- 
tron star might have passed a dense fragment of the CrA 
molecular cloud a few thousand years ago which could have 
affected a potential circumstellar disc, as well as have en- 
abled a brief history of accretion form the ISM. 
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